Depletion driven adsorption of colloidal rods onto a hard wall 
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In a mixed suspension of rods and small polymer coils, the rods adsorb onto a hard wall in 
contact with the suspension. This adsorption is studied in the low density of rods limit. It is driven 
by depletion forces and is much stronger for long rods than for spheres. This is shown by means of 
exact, numerical, calculations and an approximate theory. 



67.70.+n, 82.70.Dd, 61.25.Hq 



I. INTRODUCTION 

Most of the properties of a colloidal suspension depend 
on the shape of its constituent particles. Certainly, the 
bulk phase behaviour of colloidal particles is sensitive to 
their shape, as well as to the strength of any attractive 
interactions. Spheres behave differently from rods which 
in turn behave differently from discs. This is true in 
the presence |^,^ and in the absence ||,^ of attractive 
interactions. However, a description of the bulk phase 
behaviour does not complete the description of a suspen- 
sion. We might also like to know, for example, its be- 
haviour near the walls of its container. We study this be- 
haviour here for rodlike colloidal particles in the presence 
of nonadsorbing polymer, and compare this behaviour to 
that of spherical particles. A pure suspension of rods 
does not adsorb onto a hard wall, in the absence of sig- 
nificant van der Waals attractions; here we will assume 
that they are negligible. Indeed, the density of rods near 
the wall is below the bulk density. However, the rods in a 
mixed suspension of rods and small nonadsorbing (onto 
either the rods or the wall) polymer coils do adsorb onto 
the wall. This adsorption is driven by the increase in 
the volume available to the polymer coils when a rod is 
near a wall, the so-called depletion forces We find 

much stronger adsorption for rods than has been found 
for spheres ]lO|-p^. 

The depletion driven adsorption of spherical colloidal 
particles onto a wall has been extensively studied pO|-p[ . 
These studies were primarily experimental but compar- 
ison was made with theories for depletion forces [ p^Jl6| 
and qualitative agreement found. Both the sphere and 
the wall exclude the centres of mass of the polymer coils 
from a volume which extends up to about the radius of 
gyration of the polymer away from the surface of the 
sphere or wall. This is true so long as the radius of gyra- 
tion is not too much larger than the radius of the sphere 
|p7| . When a sphere approaches the wall, the volume 
the sphere excludes to the polymer and the volume the 
wall excludes to the polymer overlap. Thus the total vol- 
ume denied to the polymer coils goes down and so their 
translational entropy goes up. For polymer coils rather 



smaller than the spheres, free energies of adsorption of up 
to dksT have been obtained |10 1^. ks is Boltzmann's 
constant and T is the temperature. 

Rodlike colloids are not uncommon, examples are the 
tobacco mosaic virus and synthetic colloidal rods 
1^,^. The effect of adding polymer to their suspen- 
sions has been studied, at least in the bulk |P, p^ , po|| . As 
far as we are aware, no experiments comparable to those 
for spheres near a wall have been performed. However, 
Buitenhuis et al ]2C| ] report a thin layer of the nematic 
phase of the rods against the wall of their capillary. This 
is at least suggestive of some attraction between the rods 
and the capillary. The adsorption of rods does not seem 
to have been much studied theoretically. An exception is 
the work of Matsuyama et al |^ , who studied a model of 
rods with discrete orientations. However, the adsorption 
of polymers which have some degree of rigidity has been 
studied, see Ref. |22) and references therein. The results 
are consistent with those found here, as the rigidity of 
the polymer is increased, so does the adsorption. 

In the following section we first derive the exact statis- 
tical mechanical expression for the density profile of rods 
near a wall, as a function of the polymer density. We use 
this expression in two ways: The first is by evaluating 
it exactly using Monte Carlo integration and the second 
is by deriving a simple analytical approximation to it. 
Example profiles of the density of rods near the wall, 
are shown and discussed in Section III. We end with a 
conclusion. Section IV. 



II. THEORY 

We start by defining our models for the rod, the poly- 
mer coil and the wall. The interactions between polymer 
coils are neglected, i.e., a fluid of just polymer coils is 
simply an ideal gas. This leaves us with just three inter- 
actions in order to describe our system: The rod-wall, 
polymer-wall and rod-polymer interactions. The rod- 
wall interaction is that of a hard spherocylinder with a 
smooth hard wall. The spherocylinder is of length I and 
diameter d — 2rc. Thus, no point on the centreline of the 
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rod, a line of length I which runs along the centre of the 
cylindrical portion of the spherocylinder, may be within 
Tc of the wall, see Fig. 0(a). Note that, as usual the 
length I is the length of the cylindrical part of the sphero- 
cylinder; its total length is I + d. For the purposes of the 
polymer-wall and polymer-rod interactions, a polymer 
coil is considered to be a hard sphere of radius jv; this is 
the so-called Asakura-Oosawa model [ p^|jl^ , |2^ , ^ . Then 
the position of a polymer coil is entirely determined by 
the position of its centre of mass. This centre of mass 
is prevented from coming within Vp of the wall and from 
fp + fc from any point on the centreline of the rod. The 
rod excludes polymer molecules from a spherocylindrical 
volume of length I and diameter 2(rc + rp). 

Now that we have defined our model we write down 
the partition function 2 for a system of one colloidal rod 
and a fluid of polymer molecules at an activity Zp, 

S = ^ ^ / drdwdr^p exp[-/3u™(r, w)] x 

J]^ exp[-l3Upy,{Yi) - l3Urp{v,UJ,Yi)], (1) 

4=1, Afp 

where the coordinates of the centre of mass of the rod 
are denoted by r, and its orientation by uj. The coordi- 
nates of all Np polymer coils are denoted by r^p and the 
coordinates of the ith polymer coil by r j . As usual (3 is 
related to the temperature T by /3 = l/ksT, where ks 
is Boltzmann's constant. The polymer's activity Zp is re- 
lated to its chemical potential fip by Zp = exp(/?/Xj,), 
where A is the thermal volume of a polymer coil. Urw, 
Upt^, and Urp are the energies of interaction of the rod 
with the wall, a polymer molecule with the wall, and the 
rod with a polymer molecule, respectively. They are all 
hard-core interactions and so are zero unless the particles 
overlap, in which case they are infinite. As there are no 
polymer-polymer interactions, Eq. (Ill) simplifies to 



S = y drdtjcxp [— /3ur-u;(r,i^)] X 



X^TTtI / driexp[-/3up^(ri) -/3urp(r,w,ri) 



(2) 



The integrand of the part of Eq. (g) which is within 
parentheses is if the polymer sphere is within rp of ei- 
ther the wall or the surface of the spherocylinder and 1 
otherwise. Therefore, the integral within parentheses is 
equal to the volume of the system which is farther than 
rp from both the wall and the rod's surface, the free vol- 
ume for the polymer V/ (r, uj). If in addition we recognise 
the sum of Eq. (|^) as being just the expansion of an 
exponential function, we have 



The volume Vf may be written as a sum of three parts, 

Vf{Y,Uj) = Vb-Vexc+Vo{v,Uj), (4) 

where Vf, is the total volume available to a polymer in 
the absence of the rod and Vexc = T^lr^ + (4/3)7rrg, the 
volume excluded by the rod to the polymer molecules; 
Te ^ fc + Tp, the radius of the excluded volume sphero- 
cylinder around a rod. Vo is the volume of overlap of the 
volume excluded to the polymer by the rod and the vol- 
ume excluded to the polymer by the wall, see Fig. |l|(b). 
Thus if the rod is far from the wall then the regions ex- 
cluded to the polymer by the rod and by the wall do not 
overlap and Vo — but if the rod is close to the wall then 
these two excluded volumes overlap and so Wq > 0. The 
total volume available to the polymer has increased. It 
is this increase in the free volume, or to put it another 
way, the reduction in the volume excluded to the poly- 
mer, which is the driving force for adsorption of the rod 
onto the wall. The adsorption is not driven by a direct 
rod-wall interaction as is usually the case. Using Eq. (Q) 
we can rewrite Eq. (||) as 

S = exp [Zp(Vb - Vexc)] X 

drdu; exp [-(3urw (r, t^) + ZpVo (r , w)] ■ (5) 

As the wall is smooth both Urw and Vo are functions only 
of the angle between the surface normal and the axis of 
the rod, 9, and of the distance between the centre of mass 
of the rod and the nearest point of the wall, h. 

Now that we have integrated over the polymer coordi- 
nates, S is a one particle partition function and so the 
probability of finding the rod at coordinates (h, 9) is pro- 
portional to the integrand of Eq. (^) . We now consider 
an ideal gas of rods, i.e., rods at a density so low that in- 
teractions between the rods have a negligible effect. Then 
the rods are independent of each other and the density 
of rods is again proportional to the integrand of Eq. (||) . 
If the bulk density of the ideal gas of rods is pb then the 
density of rods p{h^ 9) is given by 



p{h, 9) pb exp [-(3urw{h, 9) + ZpVo{h, 9)] 



(6) 



S — J drdtj exp [— /?M„o(r, 



)+ZpVf{r,u;)]. (3) 



This is exact in the low density of rods limit and is the 
equation we will use to determine how strongly the rods 
adsorb onto the surface, in the presence of the polymer. 
We are interested in the h not the angle dependence of 
the density of rods and so we integrate Eq. (^) over 9 

1 

pW ^ Pb-^ J d{cos9)exp[-(3urw{r,9) + ZpVo{r,9)] , 

(7) 

where the factor of a half is the normalisation of the in- 
tegration over 9. Note that Eq. (|^) has exactly the same 
form for rods near a wall with a direct rod-wall attrac- 
tive interaction; the term —ZpkBTvo{v, 9) is an effective 
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rod-wall attractive interaction. It is attractive because 
if the rod is far from the wall vq — and then as the 
rod nears the wall vq becomes positive, thus tending to 
increase the density of rods. Equation ^ is directly re- 
lated to the change in free energy AF{h) when a rod is 
brought from within the bulk to a height h 

AF^~kBT\n{p{h)/pb). (8) 

This quantity has actually been measured experimentally 
for spheres ||l^,|ll|. As all the interactions are athermal 
the energy is zero and the free energy A_F is simply an 
entropy times T. Finally, from Eq. (|^) the adsorption T 
is easily calculated using 




dh{p[K) - pt). (9) 



just a random number between —1/2 and 1/2 and its co- 
ordinates perpendicular to the axis are those of a random 
point within a disc. In the endcaps the point is obtained 
by first generating a random point within a sphere. Then 
if the coordinate of this point along the axis of the sphe- 
rocylinder is positive 1/2 is added to it and if it is negative 
1/2 is subtracted. Generating a point randomly within a 
sphere or disc is done with a routine of Appendix G of 
Ref. 1^ . Now we have a randomly chosen point within 
the spherocylinder, in a frame fixed on the spherocylin- 
der. The spherocylinder is then rotated to an angle of 9 
with the normal of the wall. If the point's height above 
the wall is then less than rp then it is within that part of 
the volume excluded by the spherocylinder which is also 
excluded by the wall and hence is part of Vo- If it is above 
rp it is not. If many such points are generated then the 
fraction of points which lie below rp is an approximation 
to the fractional overlap Vo{h, 0)/vexc- 



A. Rod near a wall 

It is instructive to consider just a hard spherocylinder 
near a hard wall; the Zp = limit of Eq. . The poten- 
tial Urw is zero as long as the rod and wall do not overlap. 
For this to be true the lowest point of the centreline must 
be at least rc above the wall, as the surface of the sphero- 
cylinder is Tc away from the centreline. The lowest point 
of the centreline of the spherocylinder is always one of its 
two ends, see Fig. |l|(a). Therefore, Urw is zero provid- 
ing both ends of the centreline are at least rc above the 
wall and is infinity otherwise. Then, for Zp = 0, we can 
perform the integration over 9 of Eq. (|^) , to obtain the 
density of spherocyhnders as a function of height 

r pb2{h - rc)/l [h - rc) < 1/2 
p{h) = <^ (10) 
[pb (h- rc) > 1/2. 

The density of centres of mass of the spherocyhnders de- 
creases linearly as the wall is approached due to the re- 
stricted orientations of a rod close to a wall. 



B. Numerical evaluation of Vo 

The overlapping excluded volume Vo{h,9) is the vol- 
ume of overlap of a spherocylinder of length I and radius 
Te with a half-space. We evaluate it numerically using 
a Monte Carlo integration scheme Briefiy, first a 

randomly located point within the spherocylinder is gen- 
erated, in a coordinate frame fixed on the spherocylin- 
der. Generating a point within a spherocylinder is done 
by first selecting either the cylindrical part or the end- 
caps. The probability of selecting the cylinder (endcaps) 
is given by the fraction of the total volume of the sphe- 
rocylinder which is part of the cylinder (endcaps) . In the 
cylindrical part, the point's coordinate along the axis is 



C. Approximate theory 

Here we derive a simple analytical approximation for 
p{h). In the presence of polymer the highest density of 
rods will be near to the wall, as then Vg is largest. There- 
fore, an approximation should be accurate in this region. 
The volume Vo is quite small, a fraction of r^, if the angle 
9 is small, i.e, the rod is nearly perpendicular to the wall. 
However, Vq is large, a fraction of Ir^, if the rod is nearly 
parallel to the wall and close to it, h < rc + rp. This 
is assuming that I ^ rc. In fact as Vo appears as the 
argument of the exponential function in Eq. (Q) for the 
density, the density of rods close and nearly parallel to 
the wall increases as exp(Z). As we shall see in the next 
section this can lead to very strong adsorption. 

So, we require an approximation accurate for rods close 
to, and hence necessarily almost parallel to the wall. In 
Fig. |2| we have plotted Vo/vcxc as a function of 9 for 
a number of heights h. For small h the spherocylinder 
is restricted to a small range of angles by the hard-core 
spherocylinder-wall interaction. In addition Vo is not a 
strongly varying function of 9, especially near the wall. 
Therefore, we approximate Vo for the whole range of ac- 
cessible values of 9 by its value for 9 — 90°. Then, our 
approximate expression for p(h) is 

f exp [zpVoih, 90°)] {h - rc) < 1/2 

p{h) = { 

[pb [h- rc) > 1/2. 

(11) 

Wo(ft., 90°) is easy to calculate. It consists of two parts, 
the volume of the cylinder of radius rc at a height h which 
is below a height of rp and the volume of a sphere of the 
same radius which is below this height. So, 
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{h 



,1/2- 



Vo{h,90°) = { +f [2rl - 3rl{h - Vp) 
^{h-rpf] 



h <re + rp 

h> re + Tp. 
(12) 



Equation ([l^ ) is restricted to Tc > Vp and so requires 
generalisation if > rc- However, the Asakura-Oosawa 
model ^, p^Jl6| of the rod-polymer interaction is best for 
small polymer coils, with a radius of gyration no larger 
than that of the colloidal particle Therefore, 
we will not perform calculations with rp > Vc and so do 
not require a more general expression. 

For h > re + Tp, Vo{h, 90°) = and the approximation 
of Eq. ([TTI ) yields the same value for p{h) as in the ab- 
sence of polymer, Eq. (p^. However, for small and large 
(near tt) 0's Vq is still non-zero (see Fig. ||) and so Eq. 
( pi] ) is an underestimate. As h approaches the arc over 
which the rod may rotate tends to zero and then the ap- 
proximation of Eq. ( pi] ) becomes increasingly accurate. 
For small rp the rod must be quite close for the excluded 
volume of a near parallel rod to overlap with that of the 
wall. So, reducing rp will improve the accuracy of Eq. 
(|ll|). Also, if I is large then the arc over which the rod 
can rotate decreases and so our approximation becomes 
more accurate. In the opposite limit: I — > 0, the ap- 



proximation of Eq. 
limit. 



11^) is exact; as it is in the 







III. RESULTS AND DISCUSSION 

Now we are able to calculate the density of rods near a 
wall both exactly via the numerical technique of Section 
IIB and via the simple approximate theory of Section 
lie. We have done so for a number of different values 
of the length to diameter ratio l/d oi the rods, and the 
activity of the rods Zp. A dimensionless polymer activity 
is required; we define the reduced activity z — Zp(2rp)'^. 
As there are no polymer-polymer interactions and the 
density of rods is so low that they take up a negligible 
fraction of the volume, the density of polymer coils is 
equal to their activity Zp So, the volume fraction 

occupied by the coils is tt/Qz ~ (l/2)z. In Figs. ^ and^, 
we show how the density of the centres of mass of rods 
near the wall depends on their length and on the concen- 



is for a polymer 
lilc Fig. H is for 



tration of the polymer. Note that Fig. 
coils of radius half that of the rod wr 
polymer coils with the same radius as the rod. The free 
energy change on bringing a rod from within the bulk to 
a height h is proportional to — lnp(/i), see Eq. (||). The 
adsorption F is plotted in Fig. |^. 

Both the exact numerical integration of Eq. (0) over 6 
and the approximate Eq. (|l^) are shown in Figs. ^ and 



^. The approximation is seen to be excellent near the 
peak in p{h). The highest density of rods occurs only a 
little more than above the wall, where the overlap of 
the excluded volumes is largest and the orientations of 
a rod are very restricted. This finding provides support 
for the related approximations made in Ref. As can 
be seen in Fig. ^ Vq varies very little over the restricted 
range of orientations available to a rod close to the wall. 
This makes our approximation almost indistinguishable 
from the exact, numerical result at these heights. Unlike 
spheres, the maximum in p is not at h = rc. As h de- 
creases, the orientational freedom of the rod decreases, 
see Eq. (|lo|). and this tends to counterbalance the de- 
pletion attraction of the rod for the wall. In the h rc 
limit the orientational freedom of a rod tends to and 
so the density of rods tends to 0. The approximation is 
much less good near the minimum in p{h), which is to 
be expected as there the rods can rotate over a much 
larger angle. The approximation consistently underesti- 
mates p(h) because it underestimates Vq. As can be seen 
in Fig. g, Vo is a minimum for 9 = 90°. 

From Fig. Hit is clear that, in the presence of polymer, 
long rods, l/d = 10, adsorb much more strongly onto 
a hard wall than do spheres. In Fig. |^ the maximum 
density of spheres is only about twice that in the bulk. 
However, the density of rods reaches a maximum of over 
20 times the bulk value. For even longer rods the den- 
sity maximum increases exponentially and becomes very 
large. The larger density maxima correspond to free en- 
ergies of adhesion to the wall which are many times /c^T; 
rods much longer than l/d = 10 will adhere to the wall 
effectively irreversibly. The sensitivity of the adsorption 
to the length of the rod is also notable at small lengths. 
In Fig. H we see that, for this value of z, the density of 
short, l/d = 5, rods is actually less than that of spheres. 
As rods, but not spheres, approach a wall they lose orien- 
tational entropy, which tends to reduce the density of the 
rods. Thus, without polymer the density of rods but not 
of spheres is reduced near a wall, see Eq. (^0|). For long 
rods and high z this is more than counterbalanced by the 
greater excluded volume of rods as compared to spheres. 
As z is increased beyond the value of 0.2 the density of 
l/d — 5 rods near the wall increases more rapidly than 
that of spheres, and eventually becomes larger than that 
of spheres. 

The dependence of the adsorption on the relative radii 
of the rod and of the polymer may be seen if Figs. ^ 
and ^ are compared. Fig. ^ shows the density for rods 
with l/d ~ 10 and polymer of the same radius, rp = rc, 
while Fig. ^ shows the density of rods with the same 
l/d and the same z — 0.2 but with polymer of radius 
half that of the rod, rp — 0.5rc. Of course, although the 
volume fraction of polymer is the same in each case, the 
number density of the smaller polymer is 8 times that of 
the larger. If the polymer's radius is the same as that of 
the rod then for z = 0.2, which corresponds to a volume 
fraction of about 10%, the adsorption is very weak; the 
density near the wall is always less than that in the bulk. 
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However, the same volume fraction of polymer of half 
the radius of the rods induces quite strong adsorption, a 
maximum density over 20 times the bulk density. 

The adsorption F, calculated using Eq. (^), is plotted 
in Fig. |5|, as a function of z for a number of values of 
Ijd. The curves were calculated using the approximate 
expression for Eq. For the calculation of the 

adsorption F our approximation is highly accurate as it 
is accurate where the density is highest. This is particu- 
larly so at high z or Z, the density in the peak near the 
wall grows exponentially with increasing z or Z, see Eq. 
(pi]). Then the adsorption is dominated by this peak. It 
is obvious from Fig. g that the adsorption of long rods is 
very strong, or equivalently that only small volume frac- 
tions of polymer are required to produce a layer of rods 
near the wall with a much higher density than the bulk. 



IV. CONCLUSION 

In a mixed suspension of rods and small polymer coils, 
the rods adsorb onto a hard wall in contact with the sus- 
pension. The adsorption is driven by depletion forces: 
The entropy of the polymer coils goes up when the rods 
lie close to and parallel to the wall. From Figs. ^ and ^ 
we see that the density of rods has a narrow peak, which 
implies that a monolayer is formed. This is expected 
as the depletion forces are short ranged. The rod feels 
the wall up to //2 + Tc + 2rp away but the interaction is 
only strong when the rod approaches to within ~\- 2rp. 
Only then is there strong overlap of the excluded vol- 
umes. This is reflected in the density profile. Just at the 
surface there is a dense monolayer of rods lying parallel to 
the wall, then at larger separations the density of rods is 
actually below that of the bulk. The peak in the density 
is strongly dependent on /; longer rods adsorb strongly 
onto the wall. Again this is unsurprising, the excluded 
volume and hence the potential increase in the entropy 
of the polymer coils increases linearly with the length of 
the rod. 

If we compare rods with spheres, we see that in the 
absence of polymer the density of rods but not of spheres 
drops as the wall is approached. This is due to the de- 
crease in orientational freedom as a rod approaches a 
wall. In the presence of polymer, the larger excluded 
volume of a rod as compared with a sphere (if they have 
the same radius) causes the density of rods near the wall 
to increase more quickly than that of rods. The density 
of long rods becomes much higher near the wall than does 
the density of spheres. 

Equation (|^) for the density has exactly the same form 
for a rod attracted to the wall by a direct rod-wall inter- 
action. So, in the presence of any attraction with a range 
approximately equal to the diameter of the rod, the den- 
sity profiles will be qualitatively the same as those of our 
Figs. H and ^. This applies to, for example, van der 
Waals attractions 



So far, we have only considered an ideal gas of rods 
adsorbing onto a wall. This was reduced to a simple one- 
particle problem which, of course, did not show a phase 
transition. Although we will not perform any calcula- 
tions which account for interactions between rods at the 
wall it is, perhaps, interesting to speculate on what will 
happen near the wall. We have seen that the density 
near the wall can become much higher than that in the 
bulk, particularly if Ijd is large. The density peak is 
narrow due to the short ranged attraction and so only a 
monolayer is formed. So, we expect a dense monolayer, 
which would resemble a system of 2-dimensional rods. 
For strong adsorption the density of this monolayer will 
be much higher than in the bulk. This immediately sug- 
gests that even at low bulk densities, i.e., densities be- 
low a bulk transition to a nematic phase | p3[ | or between 
two fluid phases [ p5|j30| , the density of a monolayer just 
above the wall may be high enough for this monolayer to 
undergo a transition to a 2-dimensional nematic phase 
P, pl|j3l| . We expect the surface layer to be in the ne- 
matic phase for densities F of order 1/W and greater. 
Crystallisation of spheres at a wall at lower bulk densities 
than required for bulk crystallisation have been observed 
|]l3| and treated theoretically . 



ACKNOWLEDGEMENTS 

It is a pleasure to acknowledge a careful reading of the 
manuscript by J. Doye. I would like to thank The Royal 
Society for the award of a fellowship and the FOM insti- 
tute AMOLF for its hospitality. The work of the FOM 
Institute is part of the research program of FOM and is 
made possible by financial support from the Netherlands 
Organisation for Scientific Research (NWO). 



[1] P. van der Schoot and T. Odijk, J. Chem. Phys. 97 (1992) 
515 

[2] R. P. Sear, Phys. Rev. E 55, 5820 (1997). 

[3] D. Frenkel, in Liquids, Freezing and Glass Transitions, 
Les Houches session, edited by J. -P. Hansen, D. Levesque 
and J. Zinn- Justin (North-Holland, Amsterdam, 1991). 

[4] M. P. Allen, G. T. Evans, D. Frenkel and B. M. Mulder, 
Adv. Chem. Phys. 86, 1 (1994). 

[5] A. P. Gast, C. K. Hall and W. B. Russell, Faraday Dis- 
cuss. Chem. Soc. 76, 189 (1983). 

[6] H. N. W. Lekkerkerker, W. C. K. Poon, P. N. Pusey, A. 
Stroobants and P. B. Warren, Europhys. Lett. 20, 559 
(1992). 

[7] W. B. Russel, D. A. Saville and W. R. Schowalter, Col- 
loidal Dispersions (Cambridge University Press, Cam- 
bridge, 1989). 

[8] W. C. K. Poon and P. N. Pusey, in Proceedings of the 



5 



International School of Physics 'Enrico Fermi' Course 
CXXIX, edited by M. Baus, L. F. RuU and J. P. Ryck- [20 
aert (Kluwer, Dordrecht, 1995). 
[9] H. N. W. Lekkerkerker, P. Buining, J. Buitenhuis, G. J. 

Vroege and A. Stroobants, Proceedings of the Interna- [21 
tional School of Physics 'Enrico Fermi ' Course CXXIX, 
edited by M. Baus, L. F. Rull and J. P. Ryckaert (Kluwer, [22 
Dordrecht, 1995). 
[10] P. D. Kaplan, L. P. Fauchoux and A. J. Libchabcr, Phys. [23 

Rev. Lett. 73, 2793 (1994). 
[11] D. L. Sober and J. Y. Walz, Langmuir 11, 2352 (1995). [24 
[12] A. D. Dinsmore, A. G. Yodh and D. J. Pine, Nature 383, [25 
239 (1996). 

[13] P. D. Kaplan, J. L. Rourke, A. G. Yodh and D. J. Pine, [26 

Phys. Rev. Lett. 72, 582 (1994). 
[14] Y. N. Ohsima et al, Phys. Rev. Lett. 78, 3963 (1997). [27 
[15] S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 

(1954). 

[16] A. Vrij, Pure Appl. Chem. 48, 471 (1976). [28; 
[17] E. Eisenriegler, A. Hanke and S. Dietrich, Phys. Rev. E 

54, 1134 (1996). [29 
[18] S. Fraden, Proceedings of the International School of [30 

Physics 'Enrico Fermi' Course CXXIX, edited by M. 

Baus, L. F. Rull and J. P. Ryckaert (Kluwer, Dordrecht, [31 

1995). [32 
[19] P. A. Buining, A. P. PhiUpse and H. N. W. Lekkerkerker, 



Langmuir 10, 2106 (1994). 

J. Buitenhuis, L. N. Donselaar, P. A. Buining, A. 
Stroobants and H. N. W. Lekkerkerker, J. Colloid In- 
terface Sci. 175, 46 (1995). 

A. Matsuyarna, R. Kishimoto and T. Kato, J. Chem. 
Phys. 106, 6744 (1997). 

C. C. van der Linden, F. A. M. Leermakers and G. J. 

Fleer, Macromoloculos 29, 1172 (1996). 

G. J. Vrocgc and H. N. W. Lekkerkerker, Rep. Prog. 
Phys. 55, 1241 (1992). 

P. B. Warren, J. Phys. France I 4, 237 (1994) 237. 

H. N. W. Lekkerkerker and A. Stroobants, Nuovo Ci- 
mento 16D, 949 (1994). 

J. R. Rowlinson and B. Widom, Molecular Theory of 
Capillarity (Clarendon Press, Oxford, 1982). 
W. H. Press, S. A. Teukolsky, W. T. VetterUng, and B. 
P. Flannery, Numerical Recipes (Cambridge University 
Press, Cambridge, 2nd ed, 1992). 

M. P. Allen and D. J. Tildesley, Computer Simulations 

of Liquids (Clarendon Press, Oxford, 1987). 

R. P. Soar, J. Phys. France II 7, 877 (1997). 

P. G. Bolhuis, A. Stroobants, D. Frcnkol and H. N. W. 

Lekkerkerker, J. Chem. Phys. 107, 1551 (1997). 

Z. G. Wang, J. Phys France 51, 1431 (1990). 

W. C. K. Poon and P. B. Warren, Europhys. Lett. 28, 

513 (1994) 



6 





FIG. 1. Schematics of a hard spherocyhnder near a hard wall, (a) The solid spherocylinder shows it at the smallest angle 
possible without overlap of the spherocylinder with the wall and the dashed spherocylinder shows it the largest angle. The 
solid arc with arrows denotes the range of angles available to the spherocylinder at that height, (b) The cross-section of the 
overlap of the volume a rod excludes to the polymer with the volume the wall excludes to the polymer is shown as the shaded 
area. The solid lines are the surface of the spherocylinder and the wall, and the dashed lines are the surfaces of the volumes 
they exclude to the polymer. 
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FIG. 2. The fraction of the volume excluded by the rod which overlaps with that of the wall, as a function of 9 (in degrees) 
for several values of h/d. Each curve is labelled by the value of h/d. For all curves, l/d = 10 and = Vc- The curves are 
plotted over the range in 6 for which the hard core of the spherocylinder does not overlap with the wall. 
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FIG. 3. The ratio of the density of rods p{h) to their bulk density pb, as a function of height h above the wall. The 2 
solid curves are from Monte-Carlo integration, the 3 dashed curves are from the analytical theory. The curves are labelled with 
their value of l/d; the dashed curve labelled is for spheres. For all curves, rp = 0.5rc and z = 0.2. The inset graph is an 
enlargement of a small region of the main graph. 
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FIG. 4. The ratio of the density of rods p{h) to their bulk density pb, as a function of height h above the wall. The 2 solid 
curves are from Monte-Carlo integration, the 3 dashed curves are from the analytical theory. The curves are labelled with their 
value of z. For all curves, rp = Vc and l/d = 10. 
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FIG. 5. The adsorption F is plotted as a function of z. Each curve is labelled with its value of Ijd. The curve labelled is 
for spheres. For all curves, = 0.5rc. 
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